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ABSTRACT 

Two consecutive rotational transitions of the long cyanopolyyne HCnN, J = 39 — > 
38 and 38 37, have been detected in the cold dust cloud TMC-1 at the frequencies 
expected from recent laboratory measurements by Travers et al. (1996), and at about 
the expected intensities. The astronomical lines have a mean radial velocity of 5.8(1) km 
s~i, in good agreement with the shorter cyanopolyynes HCyN and HCgN observed in this 
very sharp-line source [5.82(5) and 5.83(5) km s~^, respectively]. The column density 
of HCiiN is calculated to be 2.8 x 10^^ cm~^. The abundance of the cyanopolyynes 
decreases smoothly with length to HCnN, the decrement from one to the next being 
about 6 for the longer carbon chains. 
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1. Introduction 

Several weak radio lines in the circumstel- 
lar shell of IRC+10216 (Beh et al. 1982) and 
the dust cloud TMC-1 (Bell & Matthews 1985) 
were originally attributed to the long carbon- 
chain molecule HCnN (cyano-deca-penta-yne), 
but a sensitive spectral survey of IRC-l-10216 
between 22.3 and 24.1 GHz (Bell unpubhshed) 
showed no evidence for this molecule at other fre- 
quencies predicted from the earlier assignment. 
This anomaly has now been explained with the 
laboratory detection of HCnN by means of a 
Fourier transform molecular-beam spectrometer 
which showed that the earlier assignment was in- 
correct (Travers et al. 1996). Although the deep 
IRC-l-10216 spectral survey also showed no ev- 
idence for spectral lines at the HCnN labora- 
tory frequencies (Bell, Feldman, & Avery 1992; 
Bell et al. 1992), the laboratory data have per- 
mitted a very deep astronomical search for this 
molecule in the narrow- line source TMC-1, and 
this has now been done with the new, very sen- 
sitive, dual-polarization GBT receiver in opera- 
tion on the NRAO0 140-foot (43 m) telescope. 
We report here the detection in TMC-1 of two 
transitions of HCnN near 13 GHz, close to the 
expected line intensity peak for a source at 10 K. 

2. Observations and Results 

The observations were done in two sessions 
with the 140-foot telescope, both during excel- 
lent weather. In 1996 December the J = 39 — > 38 
line at 13186 MHz was observed for a total of 35 
hours and in 1997 January the 38 — > 37 line at 
12848 MHz was observed for a total of 36 hours. 
The observations during the second session were 
undertaken remotely from Ottawa by means of 
the 140-foot DISPLAY software. 

Telescope pointing difficulties encountered in 

^The National Radio Astronomy Observatory (NRAO) is 
operated by Associated Universities Inc., under agreement 
with the National Science Foundation. 



both sessions sometimes required pointing checks 
every 30 minutes. For this reason, during most 
of each session a nearby transition of HCgN was 
observed simultaneously in the second correlator 
bank, which is available for independent observa- 
tions even when the dual-polarization receiver is 
used (provided that the auto correlator is set in 
the 4-bank mode). In both sessions, on starting 
observations each day the J = 5 — > 4 transition 
of HC5N at 13313 MHz and the 12 11 transi- 
tion of HC7N at 13536 MHz were also measured. 
To calculate accurate line intensities we also ob- 
served HC5N and HC7N near the meridian, so 
telescope gain changes could be estimated and 
corrections applied during data reduction. 

For the HCnN observations, three different 
local-oscillator (LO) frequencies were used over 
the two sessions so that features that moved with 
the LO frequency could be recognized. With the 
narrow line widths in TMC-1 (0.5 km s~^, or 
20 kHz at 13 GHz), spectral features of terres- 
trial origin are smeared over several line widths 
during the 4-5 day observing sessions. Sharp 
features that are present in each night's aver- 
age, or in both the first and last halves of the 
data, must then be either source-related or an 
artifact introduced by the correlator. Since the 
HCiiN spectra were in banks 2 and 4 in De- 
cember and banks 1 and 3 in January, U lines 
that occur in different channels in the two ses- 
sions might still be correlator related. Because 
of the change of banks between the two sessions, 
however, the appearance of the HCnN line at 
the expected position in both sessions is highly 
unlikely to be of instrumental origin. Although 
spectral baseline ripple is not usually a problem 
when frequency switching with small offsets (±29 
kHz) on narrow-line sources, the normal ibA/8 
focus modulation that removes ripple produced 
by residual standing waves between the primary 
and secondary reflectors was also employed as a 
precaution. 

In December the observed frequencies were 
near the center of the receiver passband and sys- 
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tern temperature was about 28 K; in January 
they were closer to the edge of the band and 
the system temperature was about 33 K. The 
telescope beamwidth (FWHM) was 2.'4 and the 
beam efficiency was 0.54 (C. Brockway, personal 
communication 1992). The observing technique 
used was small-offset frequency switching; the 
frequency-switched images were later removed 
as described previously (Bell, Avery, &; Watson 
1993; Bell 1996; Bell 1997). The autocorrclator 
was configured with 4 banks (2 for each polariza- 
tion), each with 256 channels. The frequency off- 
set used for HCgN and HCnN was ±6 channels. 
An offset of ±16 channels was used for HC5N 
and HC7N because of the hyperfine splitting of 
the HC5N transition. Taking both polarizations 
into account, the effective on-source integration 
time for HCnN was about 70 h in each session. 

Spectra containing the two observed transi- 
tions of HCiiN, and the J = 23 — > 22 transi- 
tion of HCgN observed simultaneously with the 
39 — > 38 transition of HCnN, are presented in 
Figure 1 plotted on a common velocity scale. All 
three spectra have been processed in an identi- 
cal fashion, which included overlapping, remov- 
ing the reference images, and Manning smooth- 
ing. The line parameters are presented in Ta- 
ble 1. The rest frequency obtained from labo- 
ratory work is included in column 2, the mea- 
sured LSR velocity is in column 3, the brightness 
temperature is in column 4, and the measured 
line width (FWHM) is in column 5. As can be 
seen from column 3, the velocities measured for 
HC7N and HCgN are close to the value of 5.85 km 
s~^ expected for this source (Kroto et al. 1978; 
Tolle et al. 1981). Although the velocities of the 
HCiiN lines have larger uncertainties because of 
the poorer signal to noise, the mean value ob- 
tained from the two transitions (5.84 km s~^) is 
in excellent agreement with the values found for 
the shorter cyanopolynes. The lines identified as 
HCiiN also have widths normal for TMC-1. 

It should be noted that in some cases the 
observing frequency used differed slightly from 



the rest frequency listed in column 1 of Table 
1. For HC5N the observing frequency used was 
13313.338 MHz. More importantly, for the HCgN 
J = 22 ^ 21 and 23 22 lines the observing 
frequencies were 12782.766 MHz and 13363.790 
MHz, respectively, and the differences between 
these and the rest frequencies in Table 1 were 
taken into account in determining velocities. For 
the J = 38 ^ 37 line of HCnN, the value ob- 
tained from the laboratory constants, 12848.731 
MHz, was used both as an observing frequency 
and in determining the velocity listed in Table 1. 

Three U lines are apparently present in the 
deep HCiiN spectra in Figure 1. For the rea- 
son mentioned earlier we are unable to rule out 
completely the possibility that these are spuri- 
ous features introduced by the correlator. If this 
is the case, however, they are seen only at very 
low levels. From the two HCnN spectra, the line 
density at this level (~ 2mK) is estimated to be 
~ 2 lines per MHz. It is easily shown that with 
this line density the chance of detecting a U line 
within one channel (4.88 kHz) of the expected 
HCiiN position is about 1%, and the probabil- 
ity of this happening for two transitions is thus 
very small: ^ 10~^. The probability of finding 
lines at the right frequencies with the right in- 
tensities and widths is clearly even less, so the 
identification of HCnN in space now seems quite 
secure. 

Spectra obtained for HC5N and HC7N (af- 
ter 4 and 20 minutes of integration, respec- 
tively) are shown in Figure 2 after overlapping. 
In these cases overlapping has been done but 
their frequency-switched images have not been 
removed. To determine the integrated line inten- 
sity of HC5N, the area under the three quadrupole 
hyperfine components was summed. To obtain 
the final single-component peak equivalent value, 
the total area was then divided by 0.4 km s~^ 
which was assumed to be the velocity width 
(FWHM) of the source (Broten et al. 1978). 

Column density estimates were made for all 
four observed cyanopolyyne molecules on the as- 
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sumption that (i) all lines are optically thin, (ii) 
all lines originate from a source of size 6.'0 by 
1/3 (Churchwell et al. 1978; Olano et al. 1988), 
and (iii) the excitation temperatures are identi- 
cal (i.e., there is a single rotational temperature). 
Calculations were carried out for rotational tem- 
peratures (Trot) of 6, 8, 10, and 12 K, following 
the analysis of Avery et al. (1992). Although it 
might be thought that molecules with the very 
large dipole moments (> 4 D) of the present 
cyanopolyynes would cool quickly and have ro- 
tational excitation temperatures lower than the 
kinetic temperature of the gas, this is not ex- 
pected to be the case for the low-lying transi- 
tions observed here. Since the Einstein A coef- 
ficients are proportional to u^i^^, the radiative 
cooling rates remain small even for transitions as 
high as J = 40, and the excitation temperatures 
are therefore expected to be close to the kinetic 
temperature. If so, the excitation temperature 
of ^ 11.5 K recently determined for HCgN in 
TMC-1 by Bell et al. (1997) is probably a good 
measure of the kinetic temperature of the gas in 
this source (c.f. Dickman 1975; Avery 1980). 

The dipole moments used in the excitation 
analyses are listed in Table 2 along with the col- 
umn densities calculated for Trot = 10 K. The 
dipole moments quoted by Snyder et al. (1986) 
for HC7N, HCgN and HCnN, which were based 
on crude ab-initio self-consistent field (SCF) cal- 
culations and modified by making empirical ad- 
justments and extrapolation, are probably too 
low by as much as 0.5 Debye. Recently, large- 
scale coupled cluster [CCSD(T)] calculations have 
been carried out for HC7N, HCgN, and HCnN by 
P. Botschwina (personal communication 1997). 
He finds /xo(HC7N) = 4.82(5) D, /xo(HCgN)=5.20(5) 
D, and /v,o(HCiiN)=5.47(5) D. These values are 
estimated to be accurate to within a few per- 
cent, even though HC7N to HCnN fall in the 
strongly non-linear part of the dipole-moment 
curve with respect to carbon-chain length. The 
column density found here for HCgN, 1.9 x 10^^ 
cm~^, is slightly lower than the value 3.2 x 10^^ 



cm^^ reported by Broten et al. (1978), although 
the dipole moment of this molecule turns out to 
be slightly lower (~ 7%) than their estimate of 
5.6 D. This is largely due to the higher antenna 
temperature that Broten et al. reported for the 
25 24 transition at 14535 MHz, which is al- 
most 50% higher than we would now predict. 

The calculated column densities are plotted 
for all four molecules in Figure 3 (upper), for ex- 
citation temperatures of 6, 8, 10, and 12 K. As 
the Figure shows, the abundance of HCnN does 
not vary significantly for excitation temperatures 
between 8 and 12 K. Figure 3 (upper) also shows 
how the abundance changes with the length of 
the cyanopolyyne chain, and how the decrements 
in abundance steepen as the excitation temper- 
ature increases. This can be seen more clearly 
in Figure 3 (lower) where, for Trot = 10 K, the 
decrements are slightly larger than the value of 

4 found previously by Broten et al. (1978) — 
they approach the earlier value, however, when 
lower excitation temperatures are assumed. 

The present detection of HCnN, a molecule 
with molecular weight (147 amu) nearly twice 
that of glycine, demonstrates that detectable 
quantities of fairly large organic molecules can 
form in astronomical sources under conditions 
markedly different from those on Earth. How- 
ever, if the column densities of these carbon 
chains continue to drop at the rate observed for 
the longest cyanopolyynes, the most intense lines 
of the next longer cyanopolyyne, HC13N, are un- 
likely to be stronger than ~ 1 mK. If so, detec- 
tion of this molecule in TMC-1 is likely to require 
the smaller, more efficient beam of the 100 m 
Green Bank Telescope now under construction 
by NRAO. 

We wish to thank George Liptak and the 140- 
foot telescope operators for their assistance, es- 
pecially during the remote observing session. We 
thank Ron Maddalena for providing the very use- 
ful DISPLAY software package required to carry 
out the remote observations. One of us (MBB) 
also thanks Pierre Brault for his assistance in 
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writing the LINECLEAN software and Jim Wat- 
son for helpful information and discussions on 
molecular physics. 
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Table 1 
Line Parameters in TMC-1^ 



Molecule 


Rest Frequency 


VLSR 




FWHM 


(Transition) 


(MHz) 


(km s~^) 


(mK) 


(km s^"*^) 


HC5N ( J = 5 ^ 4) 


13313.3'= 




1770(43) 


0.4<^ 


HC7N (12 ^ 11) 


13535.99P 


5.82(5) 


475(13) 


0.53(2) 


HCgN (22 21) 


12782.769^ 


5.84(5)s 


81.3(15) 


0.40(1) 


HCgN (23 22) 


13363.800^ 


5.84(5)g 


82.2(22) 


0.38(1) 


HCiiN (38 37) 


12848.728'' 


5.73(10) 


9.4(17) 


0.36(10) 


HCiiN (39 ^ 38) 


13186.853*^ 


5.96(10) 


5.4(7) 


0.56(10) 


U(13186.46) 




5.85 


4.8(7) 


0.26(10) 


U(13186.98) 




5.85 


6.3(11) 


0.26(7) 


U(12848.48) 




5.85 


7.4(20) 


0.26(6) 



aa^95Q = 04*^ 38"^ 39.3% dig^o = 25° 35' 36" 

^Tb = {Tp^/rjB (uncorrected for beam dilution) 

^Prom Alexander et al. 1976 

'^Assumed value 

^Prom Kroto et al. 1978 

^Prom Travers 1996 (unpublished) 

s After adjustment for the observing frequency used (see text) 
''From Travers et al. 1996 
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Table 2 

Cyanopolyyne Column Densities in TMC-1 For Trot = 10 K 



Molecule 


Dipole Moment (//q) 


Nl (10^^ 


Abundance 




(Debye) 


(cin-2) 


Decrement 


HC5N 


4.33*^ 


330 




HC7N 


4.82'^ 


110 


3.0 


HCgN 


5.20= 


19 


5.8 


HCiiN 


5.47= 


2.8 


6.8 



'^Adopting a beam-dilution factor of (l.'3/2.'4)= 0.54 (because the 
source is extended in one dimension; see Broten et al. 1978) 

^From Alexander et al. 1976 

=From P. Botschwina 1997 (personal communication) 
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Fig. 1. — Spectra in TMC-1 of two lines of HCnN above one of HCgN. Obtained by small-offset frequency 
switching, the spectra have been overlapped and the frequency-switched images removed. 

Fig. 2. — Spectra (upper) of the J = 5 ^ 4 transition of HC5N, showing partially resolved quadrupole 
hypcrfinc structure, and (lower) of the 12 — > 11 transition of HC7N. The frequency-switched images have 
not been removed. 

Fig. 3. — (top) Plot of the logarithm of the calculated column densities of the cyanopolyynes in TMC-1 
from HC5N to HCiiN for the indicated excitation temperatures. (bottom) Plot of the decrements between 
successive odd-numbered cyanopolyynes in TMC-1 for the indicated excitation temperatures. 
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